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WEAK  AND  STRONG  IGNITION 
II.  SENSITIVITY  OP  THE  HYDROGEN-OXYGEN  SYSTEM 
I.  Introduction 

In  a  previous  paper  [1],  refered  to  below  as  Paper  I,  we  compared 
data  from  experiments  on  reflected  shocks  in  hydrogen-oxygen-argon 
mixtures  [2]  to  corresponding  detailed  numerical  simulations.  We  found 
that  simulations  of  the  conditions  behind  the  reflected  shock  in  what 
is  called  the  strong  ignition  regime  of  the  pressure-temperature  plane 
gave  results  which  are  in  close  quantitative  agreement  to  those  ob¬ 
served  experimentally.  A  fairly  sharp,  uniform  reaction  wave  starts 
at  the  reflecting  wall  after  the  chemical  induction  time  has  elapsed. 
Simulations  for  physical  and  chemical  conditions  which  fall  in  the  weak 
ignition  regime  also  displayed  behavior  analogous  to  experiments.  Here 
the  reaction  wave  starts  at  some  time  earlier  than  the  chemical  induction 
time  at  a  location  away  from  the  reflecting  wall.  Analyses  of  the 
results  indicate  that  simulations  of  both  the  strong  and  weak  cases 
had  very  similar  numerical  perturbations.  However,  the  weak  ignition 
system  was  much  more  sensitive  to  nonuniformities  existing  behind  the 
reflected  shock. 

The  differences  in  behavior  called  "strong"  (or  "sharp")  ignition, 
and  "weak"  (or  "mild")  ignition,  were  first  noted  by  Soloukhin  and 
Strehlow  and  their  coworkers  [3, 4, 5, 6, 7, 8]  (See  primarily  Strenlow  and 
Cohen  [3]  and  Voevodsky  and  Soloukhin  [8]).  The  later  work  of  Meyer  and 
Oppenheim  [9]  showed  that  the  reaction  wave  started  at  some  location 
away  from  the  reflecting  wall  in  the  weak  ignition  case.  Both  Voevodsky 
and  Soloukhin  and  Meyer  and  Oppenheim  used  their  experiments  to  try  to 
delineate  a  universal  dividing  line  on  the  pressure-temperature  plane. 
Meyer  and  Oppenheim  [9]  and  Borisov  [10]  pointed  out  the  sensitivity  of 

certain  systems  to  fluctuations. 
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In  this  paper  we  Identify  and  isolate  those  mechanisms  which  we 
believe  cause  the  sensitivity  of  certain  systems  to  sound  wave  and 
entropy  perturbations.  First  we  show  that  chemical  sensitivity  can 
be  represented  in  terms  of  derivatives  of  the  chemical  induction  time 
with  respect  to  temperature  at  constant  entropy  and  at  constant  pressure. 

We  then  study  in  detail  the  effects  of  sound  wave  perturbations  on  both 
a  weak  and  a  strong  ignition  case.  This  study  reveals  that  in  the  weak 
ignition  case,  the  time-rate-of-c'nange  of  the  radicals  H,  0,  and  OH  show 
an  increase  in  time  which  is  faster  than  exponential  for  a  long  period  of 
time  prior  to  ignition.  Thia  ia  contrasted  to  the  classical  exponential 
behavior  exhibited  by  the  strong  ignition  case.  We  then  use  this 
information  to  derive  a  simplified  but  quantitative  induction  time 
model  which  gives  criteria  for  the  maximum  variations  produced  in  induc¬ 
tion  times  in  a  system,  At^v>  given  the  amplitude  and  frequency  of 
sound  wave  perturbations. 

We  also  observe  that  the  presence  of  a  perturbation  may  cause  ignition 

to  occur  first  at  some  location  which  is  heated  at  a  later  time  by  the 

shock  than  other  locations.  This  causes  what  appears  to  be  spotty  ignition. 

We  show  that  the  conditions  which  allow  ignition  to  be  spotty  depend  only 

on  the  state  of  the  system  before  any  energy  is  released.  A  communication 

effect  due  to  the  propagation  of  pressure  pulses  o.ice  energy  has  been 

released  would  tend  to  decrease  Ax  .  We  then  combine  the  calculations 

max 

of  ^Tmax  and  the  criterion  for  spotty  ignition  to  explain  the 
experimentally  observed  weak  and  strong  ignition  phenomena. 
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On  a  slightly  different  but  related  track,  Toong  and  coworkers 
have  done  extensive  work  showing  the  effect  of  nonequilibrium  chemical 
kinetics  on  the  propagation  of  sound  waves  [11,12,13,14,15]  Tcong  [11] 
has  shown  theoretically  that  acoustic  waves  may  e  amplified  or  attenuated 
due  to  changes  in  mean  conditions  resulting  from  chemical  reactions  and 
due  to  fluctuations  in  reaction  rates  due  to  acoustic  waves.  In  this 
paper  we  emphasize  the  effects  of  sound  wave  perturbations  on  chemical 
kinetics  and  show  how  these  perturbations  may  cause  confusion  in  the 
interpretation  of  experiments. 

II.  Sensitivity  of  the  Chemical  Induction  Time 

Previous  work  on  an  induction  parameter  model  [16]  has  required 
that  we  evaluate  the  chemical  induction  time,  T^,  for  the  ratio 
H2:0^:X/2:1:4  where  X  is  diluent.  The  curves  on  Figure  1,  which  show 
(2:1:4)  for  a  wide  range  of  pressures  and  temperatures,  were  calcu¬ 
lated  by  integrating  the  ordinary  differential  equations  describing 
the  chemical  rate  scheme  given  in  Paper  I.  This  scheme  has  been  ex¬ 
tensively  tested  against  experiments  and  shown  to  predict  irducuion 
times  accurately  [17] .  We  note  that  in  this  paper  the  induction 
time  is  defined  as  that  time  at  which  there  is  a  noticeable  (~20°K) 
rise  in  initial  temperature. 

Subsequently  we  have  found  that  the  induction  times  for  many  other 
dilutions  and  stoichiometries  (where  the  amounts  of  I^C^Ar  are  in  the 
ratios  of  a:b:X)may  be  obtained  to  excellent  approximation  from  the 
formula 
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.u ,«aJ  j 


(1) 


Ti  (a,b,X)  -  R(a,b,X)  ^(2:1:4) 


where 


.  1  fcu+b+X)  .  a+b+xl 

R(«,b,X)  -  J  J  . 


Teats  of  the  predictions  of  Equations  (1)  and  (2)  against  detailed  integration 
of  the  chemical  rate  equations  indicate  that  this  model  works  well  for 
the  three  cases  in  which  we  will  be  interested  here,  2:1:0,  2:1:7  and 
8:2:90.  The  first  case,  2:1:0,  is  the  case  for  which  Voevodsky  and  Soloukhin 
[8],  and  Meyer  and  Oppenheim  [9]  did  their  shock  tube  studies  of  weak  and 
strong  ignition.  The  second  two  cases,  2:1:7  and  8:2:90,  are  those  used  in 
the  shock  tube  experiments  and  detailed  numerical  simulations  in  Paper  I. 
Figures  2a,  b  and  c  show  contours  of  x^  on  the  pressure  temperature 
plane  in  the  range  0.5  <  P  2*  6  atm  and  900  <  1100°K,  Taken  together 

these  show  the  expected  effect  of  diluting  the  reactive  gases. 

The  data  in  Figure  2  were  then  used  to  evaluate  the  quantities 

-T  3Ti  I  ,  -T  9Ti  |  (3a) 

Ti  2T  JP  nd  Tt  3T  |S 


which  reflect  the  sensitivity  of  to  entropy  (temperature)  variations 
at  constant  pressure  P  and  sound  waves  at  constant  entropy  S.  The  constant 
pressure  derivative  was  round  directly  from  the  data.  However,  the  constant 
entropy  derivative  is  a  combination  of  two  terms: 


_  3t. 

T _ i 

T  3T 


y_ 

Y-l  T  8? 


where  y  is  the  ratio  of  specific  heats,  Cp/Cv.  The  second  term  on  the 
right  hand  side  of  Equation  (3)  is  always  negative.  Although  the  first 
term,  as  seen  in  Figure  (1),  may  be  either  positive  or  negative,  the  weighted 


sum  is  always  negative.  Figure  (3)  shows  the  constant  entropy  derivative 
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Fig.  2a  —  Contours  of  constant  chemical  induction  time  (in  seconds)  as  a  function  of 
temperature  and  pressure  for  the  mixture  Ha:Oa/2:l 


Fig.  2b  —  Same  as  Fig.  2a  for  H2:02:Ar/2:l:7 


i _ i _ I _  _J _ 1 

900  950  1000  1050  1100 

Temperature  IK) 

Fig.  2c  —  Seme  as  Fig.  2a  for  Hg:O2'‘Ar/8:2:90 
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900  950  1000  1050  1100 

Temperature  (K) 

Fig.  3  —  The  solid  lines  are  contours  of  the  sensitivity  parameter  to  perturbations 
at  constant  entropy  (sound  w  ,ves).  The  line  of  open  circles  indicates  the  extended 
second  limit.  Ihe  dashed  lir  indicates  the  division  between  weak  and  strong  igni¬ 
tion  derived  by  Meyer  and  Oppenheim  [9] . 
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for  the  2:1:0  case  and  Figure  4  shows  the  constant  pressure  derivative. 

We  note  that  high  values  of  the  derivative  indicate  that  the  system 
is  more  sensitive  for  those  values  of  pressure  and  temperature. 

The  results  for  the  2:1:7  and  8:2:90  cases  are  not  shown  because 
they  are  virtually  Identical  to  those  for  2:1:0.  Equations  (1)  and  (2) 

In  fact  predict  the  independence  of  the  sensitivity  derivatives  (3a) 
and  (3b)  to  variations  in  stoichiometry  and  dilution  and  hence  predicts 
their  universality,  This  report  has  bean  tasted  against  detailed  calcu- 
let ions  in  a  number  of  regimes  and  agreement  is  excellent.  However,  the 
universality  can  only  be  derived  under  the  assumption  that  binary  reas> 
tlons  dominate  the  determination  of  the  induction  time. 

Superimposed  on  these  contours  are  the  criteria  for  separating  weak 
and  strong  ignition  given  by  Voevodsky  and  Soloukhin  [8]  and  Meyer 
and  Oppenhaim  [?].  Voevodsky  and  Soloukhin  suggested  that  the  criterion 
is  related  to  the  extended  second  limit,  which  is  defined  by  the  com¬ 
petition  between  chain  branching  and  termolecular  recombination  of 
H  and  02<  Thus  a  line  is  defined  on  the  pressure-temperature  plane  at 
2k2  ■  k^,  where  k2  and  k^  are  the  chemical  reaction  rates  for  H  +  02  + 

OH  +  0  and  H  +  02  +  M  -*•  H02  +  M,  respectively.  The  criterion  suggested 
by  Meyer  and  Oppenheim  is  based  on  their  experiments.  We  note  that  both 
criteria  separate  more  sensitive  from  less  sensitive  regions,  although  it 

is  not  obvious  how  to  derive  these  results  -frdm  the  chemical  sensitivities 
3hown . 

III.  Detailed  Simulations  of  Sound  Wave  Perturbations 

In  order  to  evaluate  how  the  sensitivities  discussed  above  effect 
ignition,  we  have  performed  a  number  of  detailed  simulations  of  a  sound 


1050 


1100 


900  950  1000 

Temperature  (K) 

Fig.  4  —  Same  as  Fig.  3  except  contours  are  for  perturbation  at  constant  pressure 
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wave  in  a  reactive  mixture.  These  used  the  one-dimensional,  time- 
dependent  NRL  reactive  shock  model  described  extensively  in  Paper  I 
and  references  [18]  and  [16].  This  model  combines  the  Flux-Corrected 
Transport  algorithm  for  convective  transport  [19,20]  with  the  CHEMEQ 
algorithm  to  solve  the  ordinary  differential  equation  describing 
the  chemical  kinetics  [21,22].  Coupling  these  processes  has  been 
described  by  Oran  and  Boris  [23].  The  chemical  rate  scheme  used  is 
the  same  one  which  tfer anted  the  induction  times  shown  in  Figure  1  [1]. 

Consider  a  uniform  region  with  homogeneous  mixtures  of  hydrogen, 
oxygen,  and  argon  at  the  temperatures  and  pressures  shown  in  Table  1, 
These  are  the  two  cases  for  which  detailed  shock  tube  simulations  were 
described  In  Paper  1.  Then  at  the  beginning  of  the  calculation  the 
velocity  perturbation  shown  In  Figure  5  Is  imposed  on  the  system  at 
each  location  x  such  that 

v(x,t-0)  -  v  sin(~—  x)  (4) 

0  Sl. 

where  v  is  the  amplitude  and  S.  is  the  half  wavelength  of  the 
sound  wave  perturbation.  We  determine  S  by  deciding  how  many  periods 
of  the  wave  we  wish  to  occur  during  a  chemical  induction  time. 

Figure  6  and  7  show  the  .results  of  simulation  of  the  two  cases 
in  Table  I.  For  these  calculations  there  were  roughly  three  periods 
of  sound  wave  oscillation  in  an  induction  time.  This  means  that  S 

la 

is  1.17  cm  for  the  strong  ignition  case  and  15  cm  for  the  weak  case. 

4 

A  relatively  large  amplitude  perturbation,  vq  ■  1x10  cm/sec,  was 

chosen.  For  the  strong  ignition  case,  shown  in  Figure  6,  the  system 

ignites  first  at  the  right  hand  wall  at  about  85  ]js,  then  at  the 

center  and  left  hand  wall  and  at  92  Us.  Thus  we  observed  that  At  , 

max 

which  measures  the  time  difference  for  ignition  throughout  the  system, 
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Table  I 


Fig.  5  —  Velocity  perturbation  used  to  initialize  the  simulation  to  test 
the  influence  of  sound  waves  on  ignition 


Temperature  (K) 
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is  7  Msec  in  this  case.  This  is  a  very  small  delay  which  would  barely 
be  observable  in  the  schlieren  photographs,  even  at  this  relatively 
large  fluctuation  amplitude.  Note  that  without  the  sound  wave  pertur¬ 
bation  the  system  ignites  at  109  us.  Figure  7  is  an  analogous  graph 
for  the  weak  ignition  case.  Here  we  observe  that  the  right  har^  wall 
ignites  at  about  160  ys.,  an  order  of  magnitude  earlier  in  time  than 
the  unperturbed  case.  The  ignition  then  generates,  a  detonation  wave 
which  propagates  back  through  the  system  (Figure  8).  Figure  9  shows 

the  results  of  perturbing  the  weak  ignition  system  by  a  much  smaller 

2 

amplitude,  v  -  5x10  cm/sec.  Here  we  observe  that  At  *  50  usee 
r  o  max 

and  the  system  first  ignites  about  75  ys  earlier  than  if  there  had 
been  no  perturbation. 

These  simulations  show  us  that  there  are  two  effects  occurring  which 
cause  various  locations  in  the  system  to  have  different  induction  times. 
The  first  effect  is  caused  by  compression  and  the  resulting  temperature 
change  due  to  the  presence  of  the  sound  wave.  The  second  effect  is  a 
communication  effect  which  only  occurs  after  energy  has  been  released. 

The  first  effect  is  clearly  present  in  Figures  6,  7  and  9.  The  second 
effect  is  only  cbvious  in  Figure  7,  as  seen  in  the  propagating  detonation 
wave  in  Figure  8.  However,  energy  release  itself  generates  sound  waves 
which  propagate  away  from  the  ignition  center  .  Thus  the  influence  of 
one  location  igniting  may  be  felt  by  its  neighbors. 

Much  of  what  has  been  said  about  the  influence  of  sound  waves  on 
ignition  also  applies  to  entropy  waves.  These  types  of  fluctuations 
generate  sound  waves  which  are  the  major  mode  of  communication  in  systems 
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Fig.  9  —  Same  as  Fig.  7  but  with  a  weaker  perturbation 
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where  thermal  conduction  is  of  negligible  importance.  As  seen  in  Figuiu  4 
as  well  as  Figure  3,  there  is  a  large  difference  in  sensitivity  between 

the  weak  and  strong  Ignition  cases. 

IV.  "Spotty"  and  "SmooLh"  Ignition 

The  contours  shown  in  Figures  2,  3,  and  4,  the  detailed  sound  wave 
simulations  shown  in  Figures  6,  7,  8,  and  9  and  the  detailed  shock  tube 
simulations  in  Paper  I  are  convincing  evidence  that  certain  systems  are 
very  sensitive  to  entropy  mode  and  sound  v*ve  perturbations.  We  proceed 
now  to  explore  how  these  sensitivities  are  related  to  the  mechanisms 
causing  weak  and  strong  ignition. 

Consider  first  Figures  10a  and  b.  The  slope  of  the  straight  lines  in  both 
figures  represent  velocities  at  which  the  reflected  shock  moves  away 
from  the  reflecting  wall  and  chus  the  rate  of  change  of  position  at 
which  material  is  first  subjected  to  the  higher  temperature  and  pressure. 

Each  fluid  elemert  has  associated  with  it  a  time,  T^,  at  which  it  begins 
to  release  energy.  The  wavy  lines  superimposed  on  the  straight  line 
represent  the  effects  of  perturbations  in  the  gas  which  locally  alter 
the  induction  time  ft  m  Thus  we  see  in  Figure  10b  that  for  large 

enough  velocities  of  the  reflected  shock,  and  large  enough  amplitude 
fluctuations  or  small  enough  wavelength  fluctuations,  ignition  can 
first  occur  at  some  location  away  from  the  wall.  It  might,  in  fact, 
occur  at  several  places  essentially  simulataneously .  This  is  the 
effect  observed  in  the  waak  ignition  calculations  and  experiments 
described  abo%e.  We  shall  call  the  situation  illustrated  in  figure  10a 
"smooth"  ignition,  because  ignition  propagates  smoothly  from  the  reflecting 
wall  even  though  the  rate  may  show  a  small  fluctuation  about  the  mean. 


- - -  -  W*v*ul 
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Fig.  10  —  Schematic  illustration  indicating  how  ignition  in  a  system 
with  small  fluctuations  may  be  smooth  (10a)  or  spotty  (10b) 

We  shall  call  the  situation  Illustrated  in  Figure  10b  "spotty*’  ignition, 
because  ignition  can  occur  in  what  appears  to  be  a  discontinou^  or  spotty 
fashion. 

Finally  we  note  that  there  is  a  communication  effect  between  fluid 

elements  which  depends  on  the  speed  of  sound  in  the  mixture.  This  will 

tend  to  diminish  but  not  totally  destroy  spotty  ignition  as  pressure 

and  temperature  increases  due  to  energy  release  are  communicated. 

For  example,  from  Figure  6  we  see  that  At  *  15  us  and  since  the  sound 

max 

4 

speed  is  ^  6x10  cm/s,  the  material  at  the  two  end  walls  ignites 
independently.  Similar  arguments  about  communication  hold  for  the 
calculations  shown  in  Figure  9,  where  At^  *  ***  40  us.  However,  it  is 
very  clear  in  Figure  7  where  Ax  ■  150  us  that  it  is  the  ignition  at 
x»SL  that  induces  early  ignition  at  x  ■  C.  This  is  due  to  the  generation 
of  a  detonation  wave  which  propagates  back  into  the  mixture  and  ignites 
it.  Thus  we  see  that  the  idealized  sound  wave  calculations  presented 
above  had  intermingled  both  the  smooth-spotty  and  the  communication 


we  determine  that  the  condition  for  spotty  ignition  is 


(7) 


The  quantity  Ax^^  is  a  time  which  depends  on  the  strength  and  frequency 
of  the  fluctuation.  The  quantity  v  is  a  property  of  the  incident  shock. 

K 

Then  from  the  sound  speed  in  the  material  and  S^,  we  can  assess  whether 
or  not  communication  effects  may  be  important. 
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V.  Evaluation  of  At 
_ max 

In  order  to  study  the  weak-strong  ignition  phenomena  in  shuck  tubes, 
we  performed  the  detailed  simulations  of  reflecting  shocks  shown  in  Paper  I. 
Then  to  Isolate  the  phenomena  we  thought  were  causing  the  weak  ignition, 
we  investigated  '.he  interactions  of  chemical  reactions  and  sound  waves. 

Now  we  carry  this  simplification  one  step  further  by  attempting  to  further 
quantify  the  mechanism  which  leads  to  the  condition  of  Equation  (7)  and  is 
responsible  for  spotty  and  smooth  ignition  due  to  sound  wave  perturbations. 
This  requires  calculating  ATmav  in  terms  of  amplitude,  frequency,  and 
properties  of  the  shocked  fluid. 

Figure  11  shows  plots  of  the  logorithm  of  d[OH]/dt  as  a  function  of 
the  logorithm  of  [OH]  for  both  the  weak  and  strong  ignition  cases  listed 
in  Table  I.  Also  shown  are  times  at  which  specific  densities  of  [OH]  are 
reached.  Since  these  curves  have  been  derived  by  integrating  only  the 
ordinary  differential  equations  describing  the  chemical  kinetics  scheme 
given  in  Paper  I,  no  fluid  dynamic  effects  are  present.  Throughout  the 
time  period  for  which  these  curves  are  plotted,  there  has  been  no  signi¬ 
ficant  energy  release  in  the  system.  First  note  that  the  values  of  [OH] 

10  13  -3 

shown  for  the  weak  ignition  case  are  between  10  and  10  cm  where 

19  -3 

the  total  number  densities  are  about  10  cm  .  The  slope  from  relatively 

early  in  the  calculation  to  the  point  of  ignition  is  very  close  to  4/3. 

10  -3 

The  behavior  below  10  cm  is  less  easily  defined  by  a  single  exponent. 

In  fact  there  appears  to  be  a  very  fast  rise,  perhaps  a  linear  phase, 
until  it  joins  the  line  of  slope  4/3  at  101®  cm  \  Systems  which  have  the 
same  entropy  (i.c.,  different  phaooo  of  the  oamo  sound  wave)  show  the  same 
qualitative  behavior,  although  they  differ  in  ignition  time  and  at  the 
specific  density  at  which  the  slope  becomes  4/3.  For  example,  in  the  case 
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Fig.  11  —  Calculated  time-rate -cf'Change  and  OH  density  as  a  function 
of  OH  density  for  the  weak  and  strong  ignition  cases  studied 


when  T  ■  1020  K  and  P  »  3.92  atm,  which  has  the  same  entropy  as  the 

10  -3 

weak  ignition  case  in  Table  I,  the  slope  becomes  4/3  when  [OH] =3x10  cm 
and  Ti  =  666  ps.  Figure  12  shows  that  [H]  and  [0]  and  to  some  extent  [H02 ] 

are  locked  into  this  4/3  behavior  also. 


Now  consider  the  strong  ignition  line  on  Figure  11.  Here  the 
slope  is  almost  always  1,  although  it  appears  to  transition  to  the 
4/3  behavior  just  prior  to  ignition.  Again  we  note  that  isentropes 
of  these  conditions  show  the  same  qualitative  behavior.  The 
slope  1  behavior  has  been  previously  described  by,  for  example, 

Schott  and  Getzinger  [24]  and  has  been  used  hy  Brokav  [25]  to 
describe  the  hydrogen-oxygen  reaction  scheme. 

In  real  time,  the  slope  1  region  is  much  earlier  than  the  slope 
4/3  region.  In  fact,  we  know  from  studying  the  computations  in 
detail  that  the  slope  1  region  is  associated  with  domination  of  the 
mechanism  by  chain  branching  reactions  such  as  0  +  H2  -*■  OH  +  H  or 
H  +  C>2  OH  +  0,  and  the  slope  4/3  region  seems  to  be  associated  with 
the  termolecular  buildup  of  H0?,  through  the  reaction  H  +  0,  +  11  -►  IIO,  +  II 


Having  observed  the  behavior  of  the  slope  of  [OH]  for  the  weak 
ignition  case  and  its  isentropes,  we  exploit  this  very  well  defined 
behavior  to  develop  a  simplified  theory  to  evaluate  the  Ax^ay  required 
in  Equation  (7).  From  Figures  13a  and  13b,  which  cxnand  the  early  tine 
behavior  of  OH  in  Figure  11,  we  know  that  we  can  write 


dn 

dt 


+  S1n  +  S 


(8) 
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Fig.  12  —  Calculated  time-rate-of-change  of  0,  H,  and  H02  as  a  function 
of  their  density  for  the  weak  ignition  case 


where  we  have  let  n  ■  [OH].  The  first  two  terms  represent  the  initial 


behavior,  the  term  represents  an  exponential  term,  and  the  cenn 

represents  the  explosive  runaway  behavior  shown  in  Figures  11  and  12.  Wc 
consider  now  only  the  weak  Ignition  case  and  first  look  at  an  equation  of 
the  form 


dn 

dt 


S 

o 


s4/3r 


4/3 


(9) 


By  choosing  two  points,  t  and  t  ,  on  the  slope  4/3  line  on  Figure  13b 

3.  b 

and  ensuring  that  Equation  9  goes  through  them,  we  find  Then  the 

data  can  be  used  to  determine  S  ^  and  in  the  regime  below  the  4/3 
slope  according  to 


dn 

dt 


—  +  S  n 
n  1 


which  has  the  solution 

2  2  2S.t 

n  (t)  *  [n  (0)  +~-i-]  e  1 

S1 


(10) 


(ID 


We  then  match  this  equation  to  the  point  t  and  pick  values  of  n(0), 

S^,  and  which  best  reproduce  the  data.  The  ratio  S  /S^  may  be  found 
by  noting  that 


2 

n 

a 


-l) 


where  f  can  be  found  by  iterating  the  equation 
1  =  [l-f(l-£o)]e^ 


(12) 


(13) 
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Fig.  13a  —  Early  time  behavior  of  the  time-rate-of-change  of  OH  as  a  function  of 
OH  for  two  sets  of  pressures  and  temperatures  which  correspond  to  constant  en¬ 
tropy.  Also  shown  is  the  analytic  fit. 
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P0  *  3.777  *  108  dynes/cm2 
*  A  A  T„  -  1020  K 

P0  *  3.997  x  10*  dynes/cm2 
®  f>  ®  Analytic  Solution  Matching 
Points 

_J _ I _ L- 

1011  1012  1013 


[OH]  (cm-3) 

Fig.  13b  —  Same  as  Fig.  13a  except  for  later  times 


with 


£ 

0 


2,2 
n  /n 
o  a 


and 


lla  % 

n  dt 
a 


(14) 


(15) 


Typical  values  of  the  parameters  required  to  solve  Equations  (9)  and  (10) 
for  the  weak  ignition  case  are  given  in  Table  II. 

Equation  (3)  is  an  ordinary  differential  equation  and  thus  represents 
only  local  variations.  Integrating  this  equation  reproduces  the  detailed 
chemical  behavior  shown  in  Figures  13.  To  use  this  simplification  we  con¬ 
sider  an  equally  spaced  set  of  points  which  may  be  described  by  these  equa¬ 
tions.  We  then  vary  the  coefficients  at  each  location  as  given  by  the  isen- 
tropic  perturbation  resulting  from  the  velocity  given  in  Equation  (4) .  Each 
point  then  reaches  values  of  OH  signifying  ignition  at  a  different  time. 

We  use  these  time  differences  to  evaluate  At  .  Thus  we  have  a  simplified, 

infix 

inexpensive  way  to  evaluate  a  which  contains  the  varying  chemical 

sensitivities  and  is  unaffected  by  communication  effects.  We  have,  in 
effect,  a  generalized  induction  parameter  model  derived  from  the  detailed 
numerical  simulations, 

VI.  Prediction  of  Spotty  Ignition 

Figure  14  shows  At  as  a  function  of  S.  for  the  weak  ignition  case 
max  L 

2  3  3 

in  Table  I  for  amplitudes  of  5x10  ,  1x10  ,  and  2x10  cm/s.  These  calcula¬ 
tions  show  that  there  is  a  maximum  value  of  At^  occuring  at  a  wavelength 
of  2Sl  »  2-5m.  At  very  short  perturbation  wavelengths  the  frequency  of 
oscillation  is  much  faster  than  the  induction  time,  so  phase  information 
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SL  (cm) 

Fig.  14  —  Graph  of  Armtx  as  a  function  of  SL  for  three  sound  wave  perturbations 
on  the  mixture  H2:O2:Ar/8:2:90  at  1000  K  and  3.72  atm.  Also  shown  are  the  re¬ 
flected  shock  velocity  which  together  with  the  At  separates  regimes  spotty  and 
smooth  ignition. 
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about  the  sound  wave  in  various  regions  of  the  system  averages  out.  The 

remaining  RMS  quadratic  affects  cause  the  left  and  right  ends  of  the  system 

to  react  slightly  faster  than  the  center  where  the  node  of  the  sound  wave 

keeps  temperature  and  pressure  constant.  This  enhanced  reactivity  is  a 

function  of  the  perturbation  amplitude,  but  not  the  frequency.  Thus  the 

At  (S..  ,  v  )  curves  go  to  a  constant  value  at  the  left.  The  slope  of 
max  l  o 

-1,  on  the  right  for  Ax max  (S^,vo)  arises  in  an  equally  straightforward 
way.  When  the  wavelength  and  periods  are  very  long,  much  longer  than  the 
induction  time,  the  temperature  and  pressure  perturbations  across  the  system 
are  very  small  and  vary  inversely  as  the  fluctuation  period.  In  this  case  the 
induction  time  difference  scales  linearly  with  Ax  ,  and  Ax  then  scales 
in  this  limit  as  1/S^.  These  mathematical  results  are  upheld  by  the  de¬ 
tailed  calculations  which  also  show  a  weak  resonance  of  the  value  of 

Ax  whenever  the  induction  time  is  an  integral  or  half  integral  multiple 
max 

o J  the  oscillation  frequency. 

The  straight  lines  on  Figure  14  correspond  to  reflected 

:ock  velocities.  This  theory  predicts  spotty  ignition  when  a  particular 

ai  wavelength-amplitude  combination  falls  to  the  left  of  the  shock 

vel  -ity  line,  and  smooth  ignition  if  it  falls  to  the  right.  Consider  for 

3 

example  the  case  of  ST  =  15  cm.  When  v  ■  1x10  cm/s,  ignition  may  be  spotty 

L  O 

4 

for  shock  velocities  less  than  about  6x10  cm/s.  As  we  can  see 

qualitatively  from  Figure  10  and  quantitatively  from  Figure  25,  a 

higher  shock  velocity  ensures  smooth  ignition.  For  the  lower  amplitude, 

5x10  cm/s,  ignition  is  smooth  to  a  velocity  of  about  2x10  cm/ sec.  We 

note  that  the  communication  effect  which  tends  to  lower  Ax  would  have 

max 
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Che  same  qualitative  effect  as  decreasing  the  amplitude  of  the  pertur¬ 
bation. 

Even  though  the  curve  for  the  strong  ignition  case  is  not  shown  in 

Figure  14  we  may  still  use  it  to  see  why  spottiness  does  not  occur. 

Since  the  Induction  time  was  shorter  than  in  a  weak  ignition  case,  the 

value  of  corresponding  to  the  same  number  of  oscillations  in  an 

induction  time  is  much  smaller,  i.e.,  1.13  cm  compared  to  15  cm.  The 

value  of  (Sl,  vq)  is  also  much  smaller  for  equivalent  amplitudes 

2 

of  the  perturbation,  about  1  ys  for  vq  ■  5x10  cm/s.  To  get  spotty 
ignition  the  reflected  shock  velocity  would  have  to  be  very  high. 

Thus  we  see  two  reasons  why  spottiness  would  not  be  detected:  1)  we 
would  require  very  high  reflected  shock  velocities  for  spottiness  to 
occur  and  2)  we  would  have  to  have  1  ys  resolution  in  the  schlieren 
photography  to  see  spottiness. 

In  both  the  strong  and  weak  ignition  cases  studied  in  Paper  I 
and  summarized  in  Table  I,  the  velocity  perturbations  were  approximately 

3 

the  same  size.  Maximum  velocities  of  about  6x10  cm/s  near  the  re- 

2  3 

fleeting  wall  decayed  quickly  to~10  -10  cm/s  or  less.  Reflected 
shock  velocities  were  4-5x10^  cm/s.  Thus  the  strong  ignition  case 
falls  far  into  the  smooth  ignition  region  and  the  weak  ignition  case 
falls  in  the  spotty  ignition  region. 
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VII.  Conclusion 

Figures  3  and  k  suggest  that  certain  regions  in  the  pressure- 
temperature  plane  are  very  sensitive  to  entropy  and  sound  wave 
perturbations.  For  a  fixed  stoichiometry,  the  more  sensitive  regions 
tend  to  have  longer  induction  times. 

The  reason  for  this  sensitivity  lies  in  the  chemical  reactions 
and  rates  controlling  the  radical  formation  even  before  aay  signifi¬ 
cant  amount  of  energy  is  released.  When  the  reaction  H  +  0^  +  M  HO^  +  M 
becomes  significant  compared  to  H  +  0^  OH  .+  0,  the  growth  of  the  radicals 
becomes  faster  than  exponential.  We  have  seen  this  in  Figures  11-13  in 
which  the  slope  of  log  [d[0H]/dt]  versus  log  [OH]  changes  from  1  to  4/3 
at  low  enough  temperatures.  We  see,  for  example,  that  a  perturbation 
on  a  constant  adiabat  when  the  slope  is  4/3  would  create  more  of  a 
change  in  induction  time  than  an  adiabatic  perturbation  on  the  slope  1 
curve . 

For  a  system  initially  at  a  given  pressure  and  temperature, 
we  can  quantify  the  sensitivity  to  a  specific  form  of  perturbation  by 
calculating  the  range  of  induction  times  it  produces.  At _ (S  ,  v  ) .  This 

max  Li  o 

can  always  be  done  using  a  detailed  simulation.  However,  for  the  specific 
case  of  sound  waves,  a  generalized  induction  parameter  model  was  derived  above 
based  on  a  parameterization  of  the  calculated  behavior  of  the  radical  OH. 

As  shown  qualitatively  in  Figure  10,  ignition  may  be  spotty  or  smooth 
depending  on  the  relationship  of  ^Tmax  t0  the  velocity  of  the  reflected 
shock  and  the  wavelength  of  the  fluctuation.  This  result  ir  shown  quantita¬ 
tively  in  Figure  14  for  the  weak  ignition  case.  Although  the  arguments  made 
above  were  made  for  sound  waves,  they  also  hold  in  general  for  entropy  per¬ 
turbations.  That  is,  it  is  the  relation  of  the  At  characteristic  of  the 

Ula2C 


perturbation  in  the  system  to  the  velocity  of  the  reflected  shock  that 
determines  if  ignition  is  spotty.  In  some  cases  the  ^TnaX  might  be 
decreased  by  the  communication  effect  discussed  above  once  energy  begins 
to  be  released. 

Weak  and  strong  ignition  can  thus  be  analyzed  in  terms  of  what 
we  have  learned  above.  Ignition  is  weak  if  the  parameters  characterizing 
the  system  fall  onto  the  spotty  side  of  a  picture  analogous  to  Figure  14 
and  if  At  .  is  large  enough  to  be  observable.  If  Ax  or  ST  is  too 
small,  ignition  will  appear  smooth.  If  the  sytem  falls  onto  the  smooth 
side  of  Figure  14,  or  Ax  is  too  small  to  be  observable,  ignition  appears 
to  be  strong.  Thus  the  cut-off  between  weak  and  strong  ignition  depends 
equally  on  the  sensitivity  of  the  chemical  induction  time  and  on  the 
characteristics  of  the  shock  tube  which  allow  a  certain  range  of  ampli¬ 
tudes  and  frequencies  of  sound  waves  and  entropy  perturbations  to  be  both 
generated  and  measured. 

The  different  criteria  proposed  by  Voevodsky  and  Soloukhin  [8]  and 
Meyer  and  Oppenheim  [9]  may  be  interpreted  in  these  terms.  Voevodsky 
and  Soloukhin' s  determination  of  the  dividing  line  in  the  pressure-tempera¬ 
ture  plane  is  the  extended  second  limit  calculated  with  particular  values  of 
k^  and  k^.  This  is  a  purely  chemical  criterion  which  is  independent  of  the 
specific  form  of  the  driving  fluctuation.  Meyer  and  Oppenheim  determined 
experimentally  that  the  dividing  line  was  at  3x^/3T  |  *  -2psK  \  which 

occurs  for  higher  temperatures  at  low  pressures.  Their  criterion  depends 
not  only  on  chemistry  but  also  on  a  particular  type  of  fluctuation  -  in 
this  case  thermal  fluctuations  but  not  sound  waves.  We  conclude  that  the 
differences  between  the  two  criteria  proposed  can  be  explained  by  different 
levels  of  sound  and  entropy  fluctuations  in  the  two  experimental  systems. 
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We  conclude  then  that  the  analysis  presented  in  this  paper  could  be 
used  to  provide  Information  about  the  nonuni fo roi ties  existing  behind 
shock  waves.  The  possibility  also  arises  of  using  controlled  perturba¬ 
tions  to  study  and  calibrate  the  sensitivity  of  reactive  systems  and 

4/3 

of  using  the  explosive  behavior  of  the  n  regime  as  an  amplifier 
in  order  to  measure  very  low  radical  concentrations. 
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